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Abstract 

We study laser-assisted photoionization by attosecond pulses using a 
time-independent formalism based on diagrammatic many-body pertur- 
bation theory. Our aim is to provide an ab inito route to the "delays" 
for this above-threshold ionization process, which is essential for a quan- 
titative understanding of attosecond metrology. Wc present correction 
curves for characterization schemes of attosecond pulses, such as "streak- 
ing" , that account for the delayed atomic response in ionization from neon 
and argon. We also verify that photoelectron delays from many-electron 
atoms can be measured using similar schemes if, instead, the so-called 
continuum-continuum delay is subtracted. Our method is general and it 
can be extended also to more complex systems and additional correlation 
effects can be introduced systematically. 



Introduction 



The temporal aspects of photoionization from many-electron systems can now 
be explored using extreme ultraviolet (XUV) attosecond pump pulses and in- 
frared (IR) femtosecond laser probe fields. Experimental studies have evidenced 
relative delays between different ionization processes of various target systems 
including both solid-state surfaces [T and atomic gases This has stimu- 

lated a number of theoretical investigations concerning the role of the long-range 
Coulomb potential [3 El E] j many-electron screening effects [H [3 E] and electron 
localization [9^ . The delayed response of atomic systems can be decomposed as 
the sum of two terms [3 [TU] : 

TA^TW + Tec, (1) 

corresponding to the Wigner-delay of the photoelectron wave packet [11] [12] 
and the continuum-continuum delay from the laser-probe process [3] . Eq. ([T]) is 
accurate for simple single-electron systems [ini US] , but its extension to many- 
electron systems has been predicted to be difficult [2] . It is of great importance 
for attosecond metrology to understand the origin of these time delays [13 [H] , 
but also to be able to account for them in the most accurate way. Numerical 
experiments, based on the propagation of the time-dependent Schrodinger equa- 
tion, serve as an important tool to estimate the delays in photoionization, but 
inclusion of all many-electron interactions is not possible at this time [3 [51 [T7] , 
except in the restricted case of helium |18| . In these demanding numerical ex- 
periments, the fast spatial extension of the photoelectron wave packet leads to 
difficulties for the analysis of the ionization process, due to the artificial bound- 
ary of the computational box. 

In this Letter we present a method to compute complex amplitudes for this 
class of above-threshold ionization (ATI) transitions in a time-independent for- 
malism. These complex amplitudes are then used to determine the atomic 
delay in laser-assisted photoionization [13l [16] (also referred to as the "streak- 
ing delay") for electrons from the outer- most n-shell in neon and argon atoms. 
Furthermore, we gauge the validity of Eq. ^ and discuss the role of multiple 
ionization channels. Correlation effects are accounted for "all orders" of single- 
particle excitations, including the non-local exchange interaction and ground- 
state correlation. Our analysis is based on the dominant class of two-photon 
processes, where one XUV photon is absorbed from the attosecond pulse and 
one IR probe photon is exchanged. The XUV photon frequencies must then 
differ by two probe photons in order to reach the same final state, as can be 
seen in Fig.[T](a) [TH]. The photoelectron is thus probed at different intermedi- 
ate kinetic energies, in close analogy with spectral shear interferometry, yielding 
information about the spectral phase- variation of the attosecond pulse [I9j , but 
also including a characteristic response of the electronic wave packet [20] . Inter- 
estingly, these two-photon amplitudes depend rather weakly on the probe-step 
of the process in the high XUV energy range, as shown in Fig. [l] (b) for the case 
of ionization from the 3p orbital in argon. The gross features of the two-photon 
process can, therefore, be identified already in the XUV-photon matrix element. 
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Figure 1: (a) Photon picture of laser-assisted photoionization from the outer- 
most orbital of argon (3p). (b) Phase (full) and normalized probability (dash) 
calculated from the two-photon matrix elements. The thick curves correspond to 
absorption of a probe photon (A), while the thin curves correspond to emission 
of a probe photon (E). The probe photon is 1.55 eV. 



e.g. a "Cooper minima" at an XUV energy of ~ 55 eV as expected from the 
3p orbital in argon ^21,. Around this energy, the phase of the matrix element 
exhibits a non-trivial behavior as the dominant ionization rate is shifted from 
the 3p — > ed to the 3p — >■ es channel. Clearly, this complex region is ideally 
suited for a quantitative test of Eq. (IT]). 



Method 

The correlated photoelectron and ion are represented using the perturbation 
diagrams shown in Fig. [2| including linear screening for the absorption of the 
XUV photon [22] . It is convenient to rewrite the two-photon transitions in Fig. [2] 
using a perturbed wavefunction (PWF), | /9n,a ), for absorption of one XUV 
photon, il, from an initial atomic orbital, | a), with energy e^. The complex 
two-photon transition amplitude can then be expressed as a one-photon dipole 
matrix element from this intermediate PWF to the final state, | s ), through 
absorption of one probe photon, cj, 

Ms,uj.nM={s\d^\pn^a), (2) 

where global energy conservation for absorption of both photons is imposed: 
eg — Ca — il+u!. The photoelectrons are detected along the polarization direction 
of the fields, namely the quantization axis, z. Note that the dipole matrix 
element in Eq. ^ is non-trivial because it describes an ATI process, i.e. a 
transition between two wavefunctions of continuum character. 
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Figure 2: (a) Absorption of XUV (pump) and IR (probe) photons by the pho- 
toelectron. (b) Direct and (c) exchange interaction between the photoelectron 
and the remaining electrons in the core, (d)-(g) Ground state correlation inter- 
actions. The XUV photon and the IR photon are indicated by the fast wiggle 
and slow wiggle respectively, while the interaction between the photoelectron 
and the remaining hole in the atom is a dashed line. 



The PWF is setup in a Hartree-Fock (HF) basis with exterior complex scaling 
in the radial dimension 1231 : 




where the complex scaling starts at r = Rc, located far away from the atomic 
core. Typical scaling parameters are: Rq ~ 100 Bohr radii and ip = 43 degrees. 
This allows for the photoelectron to reach its asymptotic form before entering the 
exterior complex scaled region, as seen in Fig. [s] (a) . A broad range of correla- 
tion effects between the photoelectron and the ion can then be accounted for us- 
ing infinite-order, many-body perturbation theory (MBPT). We construct self- 
consistent equations for the forward-propagating (+) and backward-propagating 
(— ) PWF using the diagrams in Fig. [T] 
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Figure 3: (a) Typical outgoing photoelectron wave packet (PWF), following 
absorption of one XUV photon, with real (red) and blue (imaginary) components 
shown. The damping of the PWF for r > Rc — 100 Bohr radii is due to the 
exterior complex scaled basis set. (b) Radial integration path used to evaluate 
the dipole transition in Eq. Note that R < Rc, which implies that the 
"break-point" occurs in the unsealed region of space. 



is the zeroth-order PWF, i.e. before interaction with the IR photon and with 
no correlation. The sum over excited states {exc) in Eqs. ( 
and continuum states. The finite size of the computational box, Rb — 140 
Bohr radii, ensures that the HF basis remains discrete, while the use of exterior 
complex scaling results in complex eigenenergies for the continuum states and, 
therefore, the correct outgoing nature of the photoelectron wave packet (inside 
the unsealed region). 

Further, in Eq. Q, the terms on the second line correspond to direct and 
exchange interactions; and the terms on the third line represent the ground- 
state correlation. We can solve Eq. Q by numerical iteration including all 
interacting core orbitals, | 6 ) , of the atom, but for the XUV photon energies of 
interest in this Letter, it is sufficient to include only the outer-most n orbitals. 
Once convergence is found, the forward-propagating PWF is used to evaluate 
the two-photon matrix element in Eq. ([2]). The backward-propagating PWF is 
only required for a self-consistent treatment of ground-state correlation effects 
in Eq. Q, c.f. Ref. [21]. Also, a PWF for the process where the IR photon is 
absorbed first, i.e. time-ordered before absorption of the XUV photon, can be 
constructed in a similar fashion, but we have verified that these processes are 
weaker and can be omitted to good approximation. 

In order to evaluate Eq. Q, we apply the closed-shell HF approximation. 
Then, we remove the monopole contribution of a vacancy in the atom, say | a ) , 
from the MBPT operator, r^2 ^ form modified HF equations that include a 
"projected-hole" potential [IS]. This results in two important consequences: (I) 
the long-range interaction between the photoelectron and the ion, due to the 
monopole term, is included already in the modified HF basis; (II) the remaining 
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part of the interactions are short-ranged. This method can be extended also 
to more complex systems, such as molecules, provided that the asymptotic form 
of the PWF will be determined by this long-range monopole term. The final 
state, I s ), is not correlated in our analysis because there are no interactions 
(dashed lines) after absorption of the probe photon in Fig. |2] The exact form 
of the final state is not critical since the absolute phase-shift cancels out when 
constructing the atomic delay [131 US] , but inter-orbital interactions after the 
probe step may, none-the-less, present an interesting case for further studies. 

In the complex scaled region, Rc < r < Rb, the photoelectron wave packet 
is exponentially decaying, as observed in Fig. [s] (a). Therefore, it can not be 
directly used for the transition to the final state in Eq. but it is possible 
to solve this problem using the analytical properties of Coulomb waves. The 
ATI matrix element is evaluated by changing the path of radial integration in 
the complex plane before entering the complex scaled region, as sketched in 
Fig. [3] (b). The first part is integrated from r = to i? using a numerical 
final state and PWF. Both functions are then matched to Coulomb functions 
at this "break-point" and the next part of the integral is evaluated from R to 
R±iC using these analytical functions. The remaining integral path, sketched 
as a dashed curve, can be neglected provided that C is large enough, typically 
a few hundred Bohr radii. We have computed the integrals for several different 
break-points in order to verify this procedure. 

Results 

In Fig. |4j we present the atomic delays for photoelectrons from the outer-most 
n orbitals of (a) neon and (b) argon atoms. We find that the atomic delays 
from the outer orbitals, 2p and 3p respectively, are mostly insensitive to the 
coupling with the inner orbitals. This demonstrates the validity of the single- 
active orbital approximation for outer electrons. In contrast, the delays from 
the inner orbitals, 2s and 3s respectively, are strongly altered by the coupling 
to the outer orbitals and it is required to venture beyond a single-active orbital. 
This inter-orbital correlation leads to an induced delay of a few attoseconds in 
the entire energy range for neon. The relative delay between the two ionization 
channels is ~ 12 as for an XUV photon energy of ~ 105 eV, which is close to 
the value of 10.2 ± 1.3 as obtained using the time-dependent R-matrix approach 
by Moore et al [S]. Intriguingly, these theoretical values are still too small to 
properly explain the 21 ± 5 as measured experimentally by Schultze et al 
Our approach can be extended to include additional correlation processes |26) 
allowing for further investigation of this discrepancy. 

In the case of argon, we observe sharp delay structures close to the "Cooper 
minima" from the 3p and 3s orbitals respectively. The delay peak from the 3s 
orbital is only observed in the correlated calculation, as predicted by Kheifets in 
Ref. H] . The atomic delay difference between the 3p and 3s orbitals is ~ 78 as 
at an XUV photon energy of ~ 37 eV, which is consistent with the experimental 
value of 100 ± 50 presented in Ref. |4j. 
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Figure 4: Atomic delays from the outer-most orbitals of (a) neon and (b) argon 
atoms using a probe of 1.55 eV. The data including inter-orbital correlation 
is outlined, while the single-active orbital approximation is marked by black 
symbols. HF ionization thresholds are shown. 
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Figure 5: The "exact" continuum-continuum delay is determined for all outer- 
most n shell electrons in neon and argon as tqc = ta — t^/ using Eq. Q. 
The good agreement with the analytical curve (CC) from Ref. [12] shows that 
a meaningful separation of the delays can be made for atomic many-electron 
systems. 



Finally, in order to gauge the validity of Eq. ([T]) , we extract the Wigner delay, 
Tw, from the PWF using scattering states for emission along the polarization 
axis. This delay is then subtracted from the total atomic delay in order to obtain 
the "exact" continuum-continuum delay, tqc = ta ~ twj for each ionization 
process, see Fig. [5] An analytical curve from the asymptotic approximation |16) 
shows good agreement for all ionization processes in neon and argon. The large 
delay peaks of several hundred attoseconds in argon amount to rather small 
deviations of tcc, which are barely visible in Fig. [5] 

Conclusions 

In conclusion, we have developed a powerful approach to include correlation 
effects in laser-assisted photoionization. The delayed atomic response, ta, is 
calculated from first principle for all outer-most n-shell electrons in neon and 
argon atoms. Future analysis of attosecond pulse structures based on RABITT 
[HI, PROOF [57] and FROG-CRAB [25] should include ta in order to avoid 
errors of several hundred attoseconds close to atomic features, such as Cooper 
minima. The calculation of the continuum-continuum transition relies on the 
asymptotic form of the intermediate wave packet, but only far away from the 
atomic core, where it is known to be an outgoing Coulomb function. In this way, 
we have demonstrated the validity of Eq. ([T]) also for correlated atomic systems. 
Since the method presented here builds on well-developed many-body methods. 
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it can also be extended to other systems where the random phase approximation 
(RPA) has been apphed previously and the general form of the asymptotic 
continuum wave packet is known. As an example, it would be interesting to 
study the role of an IR probe field present in photoionization of N2, where 
calculations have been performed within the RPA covering the highest four 
occupied molecular orbitals [^ . 
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